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1.  Introduction  and  Definitions 


“Why  don’t  you  call  it  entropy... nobody  knows  what  entropy  really  is, 
so  in  a  debate  you  will  always  have  the  advantage.  ” 

John  von  Neumann 
From  John  von  Neumann  advice  to  Claude  Shannon 


From  a  purely  thermodynamic  point  of  view,  the  hydrocode  approach  means,  by  definition,  that 
the  specific  volume  V  and  the  specific  internal  energy  E  are  chosen  as  independent 
thermodynamic  parameters  of  a  two-parameter  substance  under  study.  Then,  the  specific 
entropy  function  S  =  S  (V,  E)  appears  to  be  the  corresponding  canonically  associated 
thermodynamic  potential  ( 1 ).  The  clause  “canonically  associated”  means  basically  two  things. 
First,  having  known  5  =  5  (V,  E) ,  any  other  thermodynamic  parameter  can  be  expressed  from  it 
by  means  of  differentiation  and  some  algebraic  operations  only.  In  particular,  no  other 
experimental  measurements  or  theoretical  data  are  required.  Second,  there  is  full  freedom  in 
choosing  the  function  5  =  5  (V,E) :  any  particular  choice  does  not  violate  the  First  and  Second 
Laws  of  thermodynamics.  In  other  words,  full  thermodynamic  consistency  is  guaranteed  when 
working  with  the  canonically  associated  thermodynamic  potential.  Of  course,  this 
thermodynamic  consistency  does  not  automatically  guarantee  the  practiced  applicability  of  the 
chosen  canonically  associated  thermodynamic  potential.  In  other  words,  the  appropriated  choice 
of  5  =  5  ( V,E )  should  be  made  based  on  comparison  with  relevant  observations  (2-4)  (also  see 
references  5-9  and  hundreds  of  references  therein). 

As  a  synonym,  the  thermodynamic  potential  5  =  5  (V,E)  will  also  be  called  the  complete  EOS 
for  hydrocode. 

Using  the  function  5  =  5  (V,E) ,  one  can  calculate  automatically  and  uniquely  the  heat  capacity 
Cv  =  CV  (V,  E)  of  the  corresponding  substance.  The  opposite  is  not  true:  having  known  the  heat 
capacity  Cv  =  Cv  (V,  A),  we  cannot  determine  uniquely  the  thermodynamic  potential  5  =  5 (V,  E ) . 
In  this  very  sense,  the  heat  capacity  function  Cv  =  Cv  (V ,E  )  is  one  of  the  examples  of  the 
incomplete  EOS  for  hydrocode.  Other  examples  of  the  incomplete  EOS  for  hydrocode  are 
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probably  more  known  to  the  readers.  Those  are,  for  instance,  the  pressure  equation  P  =  P(V,E) 
or  the  classical  Mie-Griineisen  EOS  for  hydrocode. 

When  using  two  or  more  incomplete  EOS  for  hydrocode,  we  can  face  different  situations.  They 
can  be  thermodynamically  consistent  with  each  other  but  remain  incomplete  even  acting  in 
concert  or  they  can  be  thermodynamically  consistent  with  each  other  but  contain  redundant 
information.  At  last,  and  this  is  the  most  common  situation,  they  can  be  thermodynamically 
incompatible.  For  instance,  the  incomplete  Mie-Griineisen  EOS  is  incompatible  with  the 
incomplete  EOS  Cv  =  Cv  (V ) ,  or  with  the  incomplete  EOS  Cv  =  Cv  (E ) ,  or  with  the  incomplete 
EOS  Cv  =Cv(T)  ,  where  T  is  the  absolute  temperature. 

Summarizing,  a  lot  of  information  about  the  complete  EOS  S  =  S  (V,E)  can  be  extracted  from 
the  function  Cv  =  Cv  (V,  E) .  Namely,  with  the  help  of  the  heat  capacity  Cv  =  Cv  (V,E),  the 
function  S  (V,E)  of  any  relevant  substance  can  be  presented  in  the  following  form: 

S(V,£)  =  A(V)+|  7 - ^ -  (1) 

£*  ]drjCvl(Tj,E)  +  A(v) 

E * 

where  A(V)  and  A  (V )  are  functions  of  specific  volume  V.  Here,  and  in  the  following,  the 
asterisk  is  used  for  thermodynamic  parameters  of  a  chosen  reference  state  (v*,E* )and  various 
thermodynamic  functions  in  the  reference  state  (in  this  sense,  one  should  interpret  the  notation 
P  or  T*). 

The  opposite  is  obviously  true.  With  any  mathematically  legitimate  choice  of  the 
functions  A(V )  and  A(V )the  formula  (1)  generates  a  thermodynamically  consistent  complete 
EOS  for  hydrocode.  We  will  call  (1)  the  heat-capacity-centered  skeleton  of  the  EOS  for 
hydrocode. 

For  each  appropriate  substance,  the  functions  A(V )  and  A  (E )  can  and  should  be  determined 
with  the  help  of  other  theoretical  means  or  by  means  of  comparison  with  additional  (to  the 
measurements  of  Cv  (V,E) )  experimental  data.  In  this  series  of  reports,  we  deliberately 
concentrate  on  the  methods  of  recovery  of  the  complete  EOS,  based  on  using  additional 
experimental  data.  To  that  end,  we  need  experimental  setups  and  appropriate  measuring 
instruments.  Sometimes,  there  are  many  available  tools,  sometimes  only  few.  Sometimes,  there 
are  no  devices  available  for  direct  measurements  (for  instance,  for  substances  in  the  interiors  of 
the  Earth  and  other  celestial  bodies).  In  the  last  situation,  the  required  information  can  be 
obtained  from  theoretical  or  indirect  measurements  (often  demanding  solution  of  inverse 
problems). 

Obviously,  for  the  recovery  of  two  functions  A(V)  and  A(E),  we  need  two  series  of 
experimental  data.  The  recovery  formulas  depend  on  the  available  experimental  setups  as  well 
as  on  the  measuring  instruments  available  or  deliberately  chosen.  Therefore,  for  one 
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substance,  there  is  a  wide  variety  of  recovery  formulas.  Next,  we  present  the  recovery  formulas 
for  some  choices  of  experiments  and  measurements.  These  formulas  give  the  readers  several 
potentially  useful  prototypes  for  direct  usage  or  as  patterns  for  the  search  the  readers’  own 
recovery  formulas. 

The  recovery  formulas  are  much  easier  when  the  heat  capacity  Cv  is  constant 

Cv  (V,  E)  =  C  =  const .  In  this  case,  both  integrals  in  equation  1  can  be  calculated  explicitly,  so  it 

reads 


S(V,E)  =  A(V)  +  Cln 


£-£, +CA(V) 
CA(V) 


(2) 


It  is  convenient  to  introduce  novel  arbitrary  functions  A(V)  and  B{  V)  instead  of  the  functions 
A(V)  and  A(V): 

A(V)  =  A(V)-C\n^p-,B(V)  =  Er  +  CT*-CA(V),  (3) 


and  to  rewrite  equation  2  in  the  form 


S(V,E)  =  A(V)  +  C\n 


E+CT*  -B(V) 
CT* 


(4) 


Using  standard  thermodynamic  identities,  we  arrive  at  the  following  formulas  of  the  absolute 
temperature  T(V,E ) : 


and  pressure  P(V,E ) 


t(v,e) 


E  +  CT*  ~B(V) 
C 


P(V,E ) 


-B'(y)-A!(v) 


E  +  CT*  ~B(V) 
C 


(5) 


(6) 


The  assumption  of  constant  heat  capacity  Cv  allows  making  many  recovery  formulas  more 
explicit  and  transparent.  It  is  both  for  deeper  penetration  into  the  corresponding  physics  and  for 
developing  benchmark  solutions  for  more  general  and  complicated  models.  However,  it  is  not 
just  a  convenient  case  for  theoretical  analysis.  Not  being  universal,  the  assumption  of  constant 
heat  capacity  Cv  is  often  successfully  used  for  gases  and  liquids.  Moreover,  experimenters 
Dulong  and  Petit  noticed  in  the  early  the  19th  century  that  this  assumption  can  be  successfully 
applied  to  solids  at  room  and  higher  temperatures.  That  is  why  we  use  the  term  the  Dulong-Petit 
(D-P)  model  of  solid  when  dealing  with  the  assumption  of  constant  heat  capacity. 
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2.  Recovery  of  the  D-P  EOS  From  Experimental  Data 


2.1  Restoration  of  the  D-P  EOS  From  the  Isentrope  Measurements 

Assume  that  we  are  capable  of  making  isentropic  loading-unloading  experiments  and  consider  a 
static  adiabat  S  =  S  passing  through  the  ambient  state.  Assume  that  we  have  devices  measuring 
the  absolute  temperature  and  pressure  in  the  range  of  interest.  Let  Ts  =  T*  (V )  and  P*  =  P*  (V ) 
be  the  values  of  the  temperature  and  pressure  at  S  =  S  .  In  terms  of  these  functions,  the  required 
operational  equation  of  state  can  be  presented  in  the  following  form: 

E-E*  +  CT*  (V)+J  dtuP •*  (or) 

Z(V,E)  =  C  in - — -f - +  5*.  (7) 

V  ’  CT*  (V) 

Using  thermodynamic  identities,  we  get  the  following  formulas  of  the  temperature  and  pressure 
functions: 

T(V,E)^T;(V)  +  ^-E  +  E\dmP;(m),  (8) 

c  cy. 

and 

/  \  *i  \  \  *  r  «/  s”|^lnJ’*(y) 

P(V,E)  =  Ps(vy  E-E  +\dmPs(m)  (9) 

2.2  Restoration  of  the  D-P  EOS  From  the  Isotherm  Measurements 

Consider  an  isotherm,  T  =  T* ,  passing  through  the  ambient  state.  Let  S*  =  S*.  ( V )  and 
P*  =  Pj  (V )  be  the  values  of  the  entropy  and  pressure  on  that  isotherm.  In  terms  of  these 
functions 

,  X  "  E-E*  S*(V)-S *  1  vt  t/  "  , 

Z(V,£)  =  Cln  1  +  — ---^U - +  —  jdmPT(m)  +ST(V).  (10) 

In  order  to  make  the  last  formula  useful  for  the  EOS  recovery,  we  need  a  calorimeter,  i.e.,  a 
device  allowing  to  measure  heat  flux.  Such  a  device  should  be  used  for  measuring  the  heat 
supply  QT  (V )  along  the  isotherm  T  =  T* .  That  measurement  would  help  to  evaluate  the  entropy 

change  S*  (V )  -  S*  =  Q*  (V )  / T*  and  rewrite  the  EOS  in  the  form  of  equation  1 1  as  follows: 
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(11) 


E-E*-Q*  (V)  +  CT  +  J  dtuP*  (m) 

*(V,E)  =  C  In - — - £ - +  S'T(V). 

Using  the  last  formula  and  thermodynamic  identities,  we  arrive  at  the  following  formulas  of  the 
temperature  and  pressure  functions: 

T(V,E)  =  T*  +  -  T*  St  ^  ~  5  +lf  dwp;  (m),  (12) 


and 


P(V,E)  =  P;(V)  + 


1  dS*(V) 
C  dV 


E-E  - 


r  [5;  (V)  -  S* ]  +  j  dmP*  (mV 


(13) 


2.3  Restoration  of  the  D-P  EOS  From  the  Hugoniot  Adiabata  Measurements 

The  presence  of  the  shock  waves  of  extremely  high  intensity  is  a  characteristic  feature  of 
numerous  applications  in  ballistics  of  penetration.  Therefore,  an  adequate  reproduction  of  the 
corresponding  shock  waves  is  of  a  paramount  importance  for  successful  modeling  efforts.  To 
that  end,  modelers  require  knowledge  of  the  EOS  for  the  parameters  that  are  close  to  the 
parameters  observed  in  the  shock-waves.  It  is  therefore  natural  to  use  the  shock-wave  related 
measurements  for  the  recovery  of  the  EOS  in  the  representative  range  of  parameters.  Next,  we 
give  the  recovery  formulas  for  some  types  of  measurements  relating  to  shock  waves. 


2.3.1  Hugoniot  Adiabata:  Internal  Energy  Density  +  Temperature  Measurements 

Assume  that  we  have  two  sorts  of  measurements  relating  to  the  Hugoniot  adiabata  passing 
through  the  reference  state  (V*,E  j .  We  begin  with  the  measurements  of  (i)  the  internal  energy 

E  =  Eh  (V )  on  the  Hugoniot  adiabata  and  (ii)  the  absolute  temperature  T  =  TH(V)  on  the 
adiabata.  The  analysis  of  experimental  measuring  of  the  function  E  =  E*H  (V )  can  be  found  in 
the  reviews  (5,  7).  Direct  measurement  of  the  absolute  temperature  T  =  Tn  ( V )  for  the  shock 
waves  of  high  intensity  is  a  difficult  engineering  problem. 

In  terms  of  these  two  functions,  E  =  E^  (V)  and  T  =  T*H  (V ) ,  the  complete  EOS  for  the  D-P 
model  takes  the  following  form: 


'  ’  CT"„(V) 


V  j 

dTU  — - r 

}■  T„(m) 


. Eh(uj)-E 
2  HK  J  t —  +  P 
uj -V 


dE*H  (cr) 
dV 


(14) 
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A  vigilance  is  required  when  the  last  formula  is  used  for  numerical  recovery  of  the  EOS.  There 
is  an  apparent  singularity  at  the  m  approaching  the  lower  limit  m  — »  V  in  the  integrand. 
Theoretically,  the  limit  is  finite  and  it  is  equal  to  P  IT  .  It  is  somewhat  more  difficult  to  handle 
this  issue  when  doing  numerical  experiments. 

2.3.2  Hugoniot  Adiabata:  Pressure  +  Temperature  Measurements 

Assume  that  the  pressure  measurements  P  =  PH  (V)  on  the  Hugoniot  adiabata  are  available  in 
addition  to  the  temperature  measurements T  =  TH(V)  and  that  the  internal  energy  measurements 
E  =  Eh(V)  become  redundant.  In  terms  of  the  two  functions,  P  =  PH(V)  and  T  =  TH(V) ,  the 
complete  EOS  for  the  D-P  model  takes  the  following  form: 


X(V,£)-S*  =  Cln 


e-e.P(p’h  +  p')(v -V)+ct’h(v) 

cr'lv) 


Pa(m)-Pl 


dpH  H 

elm 


(gt  -E*) 


(15) 


2.3.3  Hugoniot  Adiabata  +  Isotherm 

When  only  one  series  of  measurements,  E  =  E*H  (V )  or  T  =  TH  (V) ,  are  available  on  the 
Hugoniot,  then  this  series  data  should  be  combined  with  additional  data  obtained  from  some 
other  experiments. 

2.3.4  The  Hugoniot  Internal  Energy  EH(V)  and  Entropy  Sj{V)  on  Ambient  Isotherm 

Assume  that  we  get  data  relating  to  the  entropy  values  S*T  (V )  on  the  ambient  isotherm 

T  =  T*.  (16) 

Then,  the  recovery  formula  reads 


^(V,£)  =  is;(V) 


+  ln- 


V  sj 

cl  me 

'.(m)-S'T(V) 

c 

clEH  (m)  2EH(m)-2E*  +  P*(m-V*y 

1 

V 

1 

12 

_ i 

J 

V * 

Sr(r)-S* 

+  CT*e  c 


CT 


(17) 


2.3.5  The  Hugoniot  Pressure  Ph(V)  and  the  Pressure  Measurements  Pi(V)  on  Ambient 
Isotherm 

This  practical  and  important  combination  of  data  is  presented  in  the  third  report  of  this  series  (4). 
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2.3.6  Hugoniot  Adiabata  +  Isentrope 

We  consider  the  case  when,  in  addition  to  one  series  of  Hugoniot  adiabat  data,  we  take  one  more 
series  of  data  relating  to  the  isentrope. 

2.3.7  The  Hugoniot  Internal  Energy  EH(V)  and  Temperature  TS(V)  on  Ambient  Isentrope 
Measurements 

Assume  that  we  get  temperature  measurements  T*  (V )  on  the  ambient  isentrope  S  =  S* .  Then, 
the  Dulong-Petit  EOS  for  hydrocode  can  be  presented  in  the  following  form: 

Z(V,£)-S* 

~C 

,  fl  £  If,  1 

=  In  < - H —  dm  — — r 

(C7y(V)  C  J  Ts(m) 


dlnTlim)  ^  .  2EH(m)-2E 

sV  ’  EH{m) - HK  J  , - P 


dm 


m-V 


(18) 


where  we  assumed  E  =  CT  . 


2.3.8  Given  the  Hugoniot  Pressure  Ph(V)  and  Temperature  TS(V)  on  Ambient  Isentrope 
Measurements 

For  this  combination  of  data,  we  arrive  at  the  following  recovery  formula: 


E(V,E)-S * 


C 


=  ln 


1  + 


E-E'P-(p'h(V)  +  P’)(V-V')  j,  p-(m)-p--d^^(m~V) 


ct;(v) 


■  +  - 
2 


1 


crs(m) 


•  (19) 
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3.  Conclusion 


The  problem  of  finding  adequate  EOS  is  one  of  the  most  important  in  many  applications.  It  is 
also  the  central  problem  of  classical  thermodynamics.  It  cannot  be  solved  by  experiment  or 
theory,  taken  separately;  combination  of  both  approaches  is  mandatory.  Thermodynamics 
discusses  how  to  make  the  combining  non-contradictory,  i.e.,  compatible  with  the  laws  of 
thermodynamics,  and  sufficient  for  avoiding  any  arbitrariness.  We  established  the  required 
explicit  formulas  for  several  sorts  of  available  experimental  data  or  data  which  might  become 
available  in  future.  We  hope  that  some  of  these  formulas  will  be  able  to  trigger  the  imaginations 
of  experimenters. 

Needless  to  say,  many  more  combinations  of  experimental  set-ups  are  possible.  They  will 
require  other  recovery  formulas.  Hopefully,  the  presented  formulas  can  play  the  role  of  the 
illuminating  prototypes  for  the  novel  ones. 
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